^^Sb-NQR study of unconventional superconductivity in the filled skutterudite 
heavy-fermion compound PrOs4Sbi2 under high pressure up to P = 3.82 GPa 
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We report ^'^''Sb nuclear quadrupole resonance (NQR) measurements of the filled skutterudite 
heavy-fermion superconductor PrOs4Sbi2 under high pressure. The temperature dependence of 
NQR frequency and the spin-lattice relaxation rate l/Ti indicate that the crystal-electric-field split- 
ting AcEF between the ground state Fi singlet and the first excited state r4 triplet decreases with 
increasing pressure. Ac-susceptibility measurements indicate that the superconducting transition 
temperature (Tc) also decreases with increasing pressure. However, above P ~ 2 GPa, both Acef 
and Tc do not depend on external pressure up to P = 3.82 GPa. These pressure dependences of 
Acef and Tc suggest an intimate relationship between quadrupole excitations associated with the 
r4 level and unconventional superconductivity in PrOs4Sbi2. In the superconducting state, 1/Ti 
below Tc = 1.55 and 1.57 K at P = 1.91 and 2.63 GPa shows a power-law temperature variations 
and are proportional to T^ at temperatures considerably below Tc. These data can be well fitted 
by the gap model A(6l) = Ao sin 6* with Aq = 3.08 fesTc and 3.04 fceTc for P = 1.91 and 2.63 GPa, 
respectively. The results indicate there exists point nodes in the gap function. 



I. INTRODUCTION 

The filled skutterudite compound PrOs4Sbi2 is the 
first praseodymium (Pr)-based heavy-fermion supercon- 
ductor with superconducting transition temperature 
= 1.85 Kii^. The heavy-electron mass was found by 
the large electronic specific heat coefficient 7 — 310- 
750 mJ/(i4r^mol)ii^i^ and the de Haas- van Alphen effect 
measurement It is confirmed that the crystal electric 
field (CEF) ground state for a Pr'^"'" ion is a Fi singlet. 

Notably, the first excited state of the T\ triplet state 
is separated from ground state of Fi singlet by a small 
gap of AcEF/fcB ~ 10 K-^i^ii^iia Because of this smaU 
Acef, the relation between the quadrupole fluctuations 



associated with the F 



(2) 



state" 



and the occurrence of 



the unconventional superconductivity has been the focus 
of discussions . 

The superconducting gap function is important for 
understanding the mechanism of unconventional super- 
conductivity in PrOs4Sbi2. First nuclear-quadrupole- 
resonance (NQR) measurement has revealed the un- 
common nature of the superconductivity^^. The spin- 
lattice relaxation rate 1/Ti shows no coherence peak just 
below Tc, nonetheless, it follows an exponential tem- 
perature dependence below Tc. Following this, a lot 
of experimental and theoretical works have been done 
to clarify the origin of uncommon superconductivity in 
PrOs4Sbi2^^'ii^'ia'20^^^^^^^ However, the 
gap structure for superconductivity in PrOs4Sbi2 has not 
yet been determined thus far. 

Our previous NQR studies on the substitution system 
Pr(Osi_^Ru:r)4Sbi2 and PrOs4Sbi2 at P = 1.91 GPa 
strongly suggested the existence of point nodes in the su- 



perconducting gap function, since Ru doping at the Os 
site in PrOs4Sbi2 as a nonmagnetic impurity induces a 
residual density of states in the superconducting gapS^ 
and 1/Ti below Tc{P) is proportional to T^ considerably 
below Tc at P = 1.91 GPa, which is well fitted by the gap 
model A{0) = Aq sin 6* with Aq = 3.08 fceTc^^. However, 
since the pressure dependence of 1/Ti below Tc{P) is 
still not confirmed yet, further NQR studies under pres- 
sure has been deserved to determine the superconduct- 
ing gap structure and to obtain further information on 
the relation between superconducting property and the 

quadrupole/magnetic fluctuations associated by the F4 
triplet state. 

In this paper, we report on the extensive ^^"^Sb-NQR 
study of PrOs4Sbi2 under high pressure up to P = 3.82 
GPa and low temperatures. Applying pressure reduces Tc 
(Ref.2) and may also change Acep^, and therefore can 
provide further information on the symmetry of the su- 
perconducting gap as well as on the mechanism of the su- 
perconductivity. We find that Acef and Tc decrease with 
increasing pressure but become almost constant above P 
~ 2 GPa up to 3.82 GPa. These pressure dependences of 
Acef and Tc indicate an intimate relationship between 
quadrupole fluctuations and superconductivity. At P = 
1.91 GPa and 2.63 GPa, 1/Ti decreases in proportion to 
T^ at very low temperatures. This strongly evidences 
that there are point nodes in the superconducting gap 
function. 



II. EXPERIMENTAL PROCEDURES 

Single crystals of PrOs4Sbi2 were grown by the Sb- 
flux method. For NQR measurements, the coarse pow- 



2 



dered single crystals are used to allow the rf mag- 
netic field to penetrate into the sample. The hydro- 
static pressure was applied by utilizing NiCrAl/BeCu 
piston-cylinder type cell for P < 3 GPa and indenter 
type clamping cell'^- for P > 3 GPa, filled with Si oil 
and Daphne 7373 as a pressure-transmitting medium, 
respectivelj«22ii22. The pressure at low temperatures was 
determined from the pressure dependence of the val- 
ues of Sn and/or Pb metals measured by a conventional 
four-terminal method. Temperature dependence of ac- 
susceptibility is measured by a four terminal method us- 
ing NQR coil. Data below 1.4 K were collected using a 
'^He/^'He dilution refrigerator at P = and a '^He refrig- 
erator under P = 1.91 and 2.63 GPa, respectively. 

The Sb nuclei have two isotopes of ^^^Sb and ^^^Sb 
with natural abundances of 57.3% and 42.7%, respec- 
tively. Since ^^^Sb and ^^'^Sb have the nuclear spin / = 
5/2 and 7/2, respectively, five Sb-NQR transitions are 
observed. In the present experiment, all measurements 
were done at the ±3/2^ ±5/2 transition (hereafter, 2vq 
transition for short) of the ^^^Sb nucleus. 

Figures. 1(a) and (b) show a typical data set of "'^^^Sb- 
NQR recovery curves of 2z/q transition to obtain relax- 
ation time Ti which are measured at T = 1.54 K just 
below Tc(P) = 1.57 K (Fig.l(a)) and T = 0.59 K weU 
below Tc(P) (Fig.l(b)) at P = 2.63 GPa, respectively. 
Since there is an asymmetry parameter 77 at the posi- 
tion of Sb nucleus, Ti is determined by the theoretical 
curve of nuclear magnetization where the value of 77 is 
incorporated^. As seen in Fig.l (b), to avoid possible 
heating due to the rf pulses, we used small amplitude rf 
pulses in the Ti measurements at very low temperatures. 
We confirmed the lack of a heating effect to Ti process 
by ensuring that the spin echo intensity is not affected by 
an rf pulse with a slightly off-resonance frequency, which 
was applied before the 7r/2 — tt pulse sequence. 



III. UNCOMMON SUPERCONDUCTIVITY IN 
PROS4SB12 AT AMBIENT PRESSURE 

Figure 2 shows the temperature dependences of l/Ti 
at ambient pressure (solid circles) . The ambient-pressure 
data are in excellent agreement with those reported pre- 
viously (open circles)^^. At high temperatures, 100 K 
> T > To ~ 10 K, 1/Ti increase slightly and become 
temperature independent with decreasing temperature. 
In this case, the relaxation rate is dominated by the 
Pr-4/^-derived localized magnetic moments on the F^'^'' 
triplet state. With decreasing temperature below Tq ^ 10 
K, l/Ti starts to decrease because a depopulation of 4/ 

electron at r4 ' triplet state occurs below T < AcEp/fcB- 
Such a phenomenon is also confirmed in temperature de- 
pendence of 2vq resonance frequency as discussed later. 
Remarkably, just below Tc, the Hebel-Slichter (coher- 
ence) peak is absent, indicating a non s-wave supercon- 
ductivity. At further low temperatures 1/Ti decreases 
exponentially.^'* Below T = 0.4 K, l/Ti become a con- 
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FIG. 1: (colour online) Recovery curves of the nuclear mag- 
netization measured at 2i^Q transition of ^^^Sb in PrOs4Sbi2 
at (a) r = 1.54 K (just below Tc{P) = 1.57 K) and (b) T = 
0.59 K (considerably below rc(P)) at P = 2.63 GPa. The red 
curves are theoretical fitting curve, (see text) 



stant, which is likely due to impurities. The present sam- 
ple reproduces the previous results, but shows smaller 
contribution of l/Ti = const, behavior at the lowest tem- 
perature region than previous report^**, which suggests 
that now samples have better quality^l. The most likely 
cause for such T- independent l/Ti is the presence of a 
small amount of magnetic impurity, which is pressure- 
independent. 

Figure 2 inset shows the temperature dependence of 
the 2uq transition at P = 0. Tq is the temperature at 
which the 2vq resonance frequency increases abruptly. 
Since the electrical field gradient (EFG) is predominantly 
determined by the on-site charge distribution, the NQR 
frequency is a powerful probe of the population of the 
ground/excited state. Indeed, in both PrOs4Sbi?i^ and 
PrRu4Sbi2^, To is in good agreement with Acef/^b- 
More recently, it has been suggested that the tempera- 
ture dependence of NQR frequency can be accounted for 
by the EFG associated with the hexadecapole moment 
of the r4^'' stated. Therefore, it is concluded that the 

increase in the NQR frequency below Tq is due to the 

(2) 

depopulation of / electron in the r4 ' state below this 
temperature. Here, we determined Tq = 12 K at am- 
bient pressure, below which l/Ti and diyq/dT starts to 
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decrease. 
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FIG. 2: (colour online) Temperature dependence of l/Ti at P 
= (solid circles) along with the data cited from literaturei^ 
(open circles). Dotted lines indicate a relation of 1/Ti = 
constant. Solid and dotted arrows indicate Tc and To, respec- 
tively. Inset shows temperature dependence of NQR reso- 
nance frequency for 2i/q transition at ambient pressure (solid 
circle) along with previous report (open circle) cited from 
literaturei^, respectively. 



IV. PRESSURE DEPENDENCE OF Acef 
PROS4SB12 



IN 



As discussed in previous section, we can estimate a 
value of Acef from the temperature dependence of NQR 
frequency. Figure 3 shows typical data set of tempera- 
ture dependence of ""^^^Sb-NQR spectra measured at 2z/q 
transition at P = 2.63 GPa. As seen in Fig. 3, aX P = 
2.63 GPa, ^^^Sb-NQR spectrum has a symmetrical shape 

and the center of spectrum clearly shifts to higher fre- 

(2) 

quency due to the depopulation of / electron in the r4 
state below Tq (cx Acef/^b)- Thus, we can determine 
the peak position of NQR spectra precisely even under 
high pressure. Figure 4 summarizes the temperature de- 
pendence of 2i'Q resonance frequency under various pres- 
sures. The temperature dependence of 2vq resonance 
frequency shifts to lower temperature under pressure in- 
dicating the decrease of Acef with increasing pressure 
but it does not change between P = 2.63 and 3.82 GPa. 
As seen in Fig. 4 inset, we define Tq at which df /dT starts 
to decrease steeply. Here, / denotes 2i'q resonance fre- 
quency. We determined Tq{P) ~ 10, 7, and 7 K at P 
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FIG. 3: (colour online) Temperature dependence of 2i^q- 
Sb-NQR spectra for PrOs4Sbi2 at P = 2.63 GPa. Solid 
curves are results of Lorentzian fittings. Solid arrows indicate 
peak positions. 



— 0.84, 2.63, and 3.82 GPa, respectively. The pressure 
dependence of Tq is summarized in Fig. 6 (a). Notably, it 
is found that Tq decreases with increasing pressure up to 
P ~ 2 GPa^^ but becomes pressure independent above 
P - 2 GPa. 
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FIG. 4: (colour online) Temperature dependence of resonance 
frequency of 2i/q transition at P = 0.84 GPa (solid circles), 
2.63 GPa (solid triangles), and 3.82 GPa (solid squares), re- 
spectively. Inset shows temperature dependence of df /dT at 
P = 2.63 GPa. / denotes 2uq resonance frequency. Dotted 
lines are eye guide. Solid arrows point to Tq. (see text) 



PRESSURE DEPENDENCE OF T 
PROS4SB12 



IN 



Next, we turn to pressure dependence of Tc in 
PrOs4Sbi2. Figure 5 shows temperature dependence of 
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ac-susceptibility at P = 0, 0.84, 1.91, and 2.63 GPa (main 
panel) measured using NQR coil in piston-cylinder type 
pressure cell and P = 3.82 GPa (inset) measured in in- 
denter type cell, respectively. As seen in the figure, the 
onset of superconducting diamagnetisms is clearly ob- 
served at all pressures. From the onset of the super- 
conducting diamagnetism, we determined Tc{P) — 1.87, 
1.73, 1.55, 1.57, 1.57, and 1.48 K at P = 0, 0,84, 1.91, 
2.34, 2.63, and 3.82 GPa, respectively. The pressure 
dependence of Tc is summarized in Fig. 6 (b). Tc in 
PrOs4Sbi2 decreases with increasing pressure^. How- 
ever, it saturates above P — 2 GPa. Notably, this pres- 
sure dependence of Tc is quite similar to the pressure de- 
pendence of To (see Fig. 6 (a)), which suggests an intimate 
relationship between Acef and the superconductivity in 
PrOs4Sbi2. 

On the theoretical side, it has been proposed that the 
superconductivity is mediated by the excitons due to the 
r4^'' — Fi quasi-quarteti^. In such case, Tc{P) would 
increase due to the enhancement of such interaction when 
Acef is reduced. However, clearly, our results indicate 
Tc decreases with increasing pressure and do not lend a 
straightforward support to the theory. Further theories 
are called for. 
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FIG. 5: (colour online) Temperature dependence of ac- 
susceptibility at P = 0, 0.84, 1.91, and 2.63 GPa, respectively. 
Inset shows an onset of superconductivity at P = 3.82 GPa. 
Arrows indicate Tc(P). 



VI. EVOLUTION OF NORMAL STATE 
ELECTRIC PROPERTIES IN PROS4SB12 UNDER 
PRESSURE 

Next, we discuss the pressure and temperature depen- 
dences of ^^■^Sb-NQR relaxation rate l/Ti, which gives 
information on the dynamical electric and on supercon- 
ducting properties. Figure 7 shows the temperature de- 
pendence of l/TiT at P = 0, 1.91, 2.34, and 2.63 GPa, 
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FIG. 6: Pressure-temperature phase diagrams for PrOs4Sbi2 
for To (solid triangles in (a)) and (solid circles in (b)) de- 
termined by present results. Dotted curves are eye-guides. 



respectively. In section IV, we concluded that Acef de- 
creases with increasing pressure and it saturates above 2 
GPa. Such a conclusion is also supported by the tempera- 
ture dependence of 1/TiT under pressure. The pressure 
effect appears below 4 K. At P — 0, the reduction of 

l/TiT results in a peak structure in the plot of l/TiT 

(2) 

versus T, which is due to the depopulation of the T\ 
state below T(^. At high pressures, 1/TiT continues 
to increase and the decreases of l/TiT occur at lower 
temperature, indicating the decrease of Acef- Since 
Acef decreases with increasing pressure, this may in- 
duce the increase of quadrupole fluctuations due to the 
f electron at the F4 state. These results are consis- 
tent with the conclusion inferred from the magnetization 
measurement^" At P = 2.34 and 2.63 GPa, l/TiT in the 
whole temperature regions has almost the same value. 
Thus, the peak position and Tc found in the temperature 
dependence of l/TiT are unchanged above P ~ 2 GPa, 
and these are consistent with the pressure dependence of 
To and Tc determined by the temperature dependences of 
2h'Q transition and ac-susceptibility as discussed in pre- 
vious sections. 
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FIG. 7: (colour online) Temperature dependences of ^^^Sb- 
NQR 1/TiT for PrOs4Sbi2 at P = (solid circles), P = 1.91 
GPa (solid triangles), P = 2.34 GPa (solid squares), and P 
— 2.63 GPa (solid diamonds) along with the data at ambient 
pressure cited from ref.M (open circles). Solid arrows indicate 
Tc{P). 



VII. SUPERCONDUCTING PROPERTIES 
UNDER PRESSURE 

Next, we discuss the superconducting property in 
PrOs4Sbi2 under pressure. Figure 8 shows the tempera- 
ture dependence of l/Ti at P = 0, 1.91, and 2.63 GPa, 
respectively. As discussed in section III, the 1/Ti is T- 
independent above Tg, indicating that the relaxation in 
the high temperature region is dominated by the Pr-4/^- 
derived localized magnetic moments,^ With decreasing 
temperature below Tq, l/Ti starts to decrease. At P = 
1.91 and 2.63 GPa, the reduction of l/Ti below Tq to 
becomes smaller compared to ambient pressure data since 
the contribution of the excitation gap Acef in 1 /Ti be- 
comes smaller under pressure. Notably, just below Tc{P), 
no coherence peak is observed at high pressures, even 
though Tc differs in each pressures. This indicates that 
the lack of coherence peak in PrOs4Sbi2 is intrinsic prop- 
erty. 

The temperature dependences of 1 /Ti below Tc (P) at 
high pressures are markedly different from that at ambi- 
ent pressure. As seen in solid line in Fig. 8, l/Ti at P 
= 1.91 and 2.63 GPa decreases in a power-law like tem- 
perature dependence below Tc{P). However, at ambient 
pressure, 1/Ti decreases exponentially as indicated by 
dotted line in Fig. 8. What is the origin of these different 
temperature dependence of 1 /Ti below Tc at ambient and 
under pressures? One possible scenario is provided be- 
low. At ambient pressure, the larger excitation gap Acef 
contributes to the reduction of l/Ti below Tc, which 
make the temperature dependence of 1/Ti exponential 
like across Tc. Actually, at P = 1.91, 2.34, and 2.63 GPa, 
the anomaly at Tc{P) in temperature dependence of l/Ti 
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FIG. 8: (colour online) Temperature dependence of Sb-NQR 
l/Ti at P = (solid circles), P = 1.91 GPa (solid triangles), 
and P — 2.63 GPa (solid squares) along with the data at 
ambient pressure cited from ref." (open circles) . The straight 
line is a guide to the eye. The dotted curve depicts the relation 
l/Ti oc exp(-Ao/fcBrc) with Ao/feTc = 3.45, proposed by 
Yogi et ali^ Arrows indicate Tc- 

becomes clearer compared to ambient pressure data since 
the Acef decreases under pressure and its contribution 
to l/Ti under pressure becomes smaller. So, we can de- 
termine Tc (P) from the temperature dependence of 1 /Ti 
as seen in Figs. 8 and 9. Such a scenario is also supported 
by our previous NQR studies in the substitution systems 
Pr(Osi_a:Rua;)4Sbi2i^ It was confirmed that Acef be- 
comes larger when Ru is substituted for Os site^^'^*. As 
a result, in Pr(Osi_j,Rua;)4Sbi2, it becomes more diffi- 
cult to determine Tc in the temperature dependence of 
1/Ti compared to PrOs4Sbi2 due to the increase of the 
contribution of excitation gap Acef in across Tc^ 
Thus, it is suggested that the excitation gap Acef /^b 
> 10 K may make temperature dependence of l/Pi look 
unusual. 

Notably, at P = 1.91 and 2.63 GPa, 1/Ti become pro- 
portional to T^ below T ~ 0.6 K. We find that a point- 
nodes model, with a low-energy(P) superconducting den- 
sity of states (DOS) proportional to E^, can explain well 
these temperature dependences of l/Pi below Tc(P). As 
plotted in Figs. 9 (a) and (b) by solid curves, we ob- 
tained excellent agreement with experimental results be- 
low Tc(P). Here, we employed Anderson-Brinkmann- 
Morel (ABM) p-wave model^i^ to represent tempera- 
ture dependence of 1 /Ti below Tc (P) as 
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where Ns{E)/No = E/VE^ - with A{0) = Aosin6'. 
From fitting results, we obtained Ao/fceTc = 3.08 and 
3.04 for P — 1.91 and 2.63 GPa, respectively. Since the 
superconducting transition temperatures at P = 1.91 and 
2.63 GPa are almost same, the values of Ao/^bTc are 
very close in both pressures. Recently, the p + h mode& 
proposed to explain the results of thermal conductivity 
would give a T'^-like dependence since the DOS at low-£' 
is linear in E, and is therefore not compatible with our 
results of 1/Ti cx T^. 
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FIG. 9: (colour online) (a) Temperature dependence of l/Ti 
at P = 1.91 GPa (solid circles in (a)) and P = 2.63 GPa 
(solid squares in (b)). Solid and dotted arrows indicate Tc 
and To, respectively. The solid curves are fit assuming the 
ABM state, (see text) The solid and dotted lines indicate the 
relation of l/Ti oc constant at high temperature region and 
1/Ti oc well below Tc, respectively. 



VIII. CONCLUSION 

In this paper, we have presented extensive ^^^Sb-NQR 
studies on the normal and superconducting states in filled 



skutterudite heavy-fermion superconductor PrOs4Sbi2 
under high pressure up to 3.82 GPa. 

From the temperature dependences of ^^'^Sb-NQR fre- 
quency and the spin- lattice relaxation rate 1/Ti under 
pressure, it is indicated that the crystal-electric-field 
splitting AcEF between the ground state Pi singlet and 

the first excited state r4 triplet decreases with increas- 
ing pressure. However, further application of pressure 
above P — i GPa reveals for the first time that Acef 
does not change in 2 GPa < P < 4 GPa. Concomitantly, 
similar pressure dependence is confirmed in the supercon- 
ducting transition temperature Tc (J^) ■ Although the ori- 
gin of these pressure dependences of Acef and Tc is still 
unknown, it suggests superconductivity in PrOs4Sbi2 is 
closely related to the r4 state. In the superconducting 
state a.t P — 1.91 and 2.63 GPa, 1/Ti shows no coher- 
ence peak just below Tc{P) and shows a power-law tem- 
perature variations and are proportional to T^ at tem- 
peratures considerably below Tc(P). These data can be 
well represented by the gap model A{6) ~ Aosin^ with 
Ao = 3.08 /cbTc and 3.04 fePc for P = 1.91 and 2.63 
GPa, respectively. This evidences superconductivity in 
PrOs4Sbi2 has point nodes in the gap function. 

These results indicate quadrupole/magnetic fluctua- 

(2) 

tions induced by P^ triplet state may play a vital role 
for the occurrence of unconventional superconductivity in 
PrOs4Sbi2. We believe present results shed further light 
on the mechanism of unconventional superconductivity 
in this compound, and will stimulate further theoretical 
works. 
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